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(*S, • Context. The CoRoT satellite has recently discovered a hot Jupiter that transits across the disc of a F9 main-sequence star called 

CoRoT-6 with a period of 8.886 days. 

Aims. We model the photospheric activity of the star and use the maps of the active regions to study stellar differential rotation and 
<""j , the star-planet interaction. 

Methods. We apply a maximum entropy spot model to fit the optical modulation as observed by CoRoT during a uninterrupted 

interval of ~ 140 days. Photospheric active regions are assumed to consist of spots and faculae in a fixed proportion with solar-like 

contrasts. 

Results. Individual active regions have lifetimes up to 30-40 days. Most of them form and decay within five active longitudes whose 

different migration rates are attributed to the stellar differential rotation for which a lower limit of AQ/Q. = 0. 12 ± 0.02 is obtained. 

Several active regions show a maximum of activity at a longitude lagging the subplanetary point by ~ 200° with the probability of a 

chance occurrence being smaller than 1 percent. 

Conclusions. Our spot modelling indicates that the photospheric activity of CoRoT-6 could be partially modulated by some kind of 

star-planet magnetic interaction, while an interaction related to tides is highly unlikely because of the weakness of the tidal force. 

Key words, stars: magnetic fields - stars: late-type - stars: activity - stars: rotation - planetary systems - stars: individual (CoRoT-6) 
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1. Introduction gions allow us to estimate a lower limit for the amplitude of 

>— » the differential rotation of CoRoT-6. We also use the maps to 

^ CoRoT (Convection, Rotation and Transits) is a photometric look for signs of star . planet interaction and find some statisti- 

space experiment devoted to asterosetsmology and the search for cal ev idence of active regions lagging the subplanetary point by 

extrasolar planets by the method of transits (Ba glin et al. || 2006|) , „ 17Q o _ 2QQO during some time intervals . The case for a star . 

It has recently discovered a hot Jupiter with an orbital period of p i anet int eraction in other systems has been recently reviewed 

8.886 days and a mass of 2.96 ± 0.34 Mj orbiting a F9 main- b y |^i |2008) and IShkolniket "^(120091 
sequence star dubbed CoRoT-6 (IFridlund et al.11201 0). The star 
3 . is magnetically active and shows a rotational modulation of its 

optical flux with a relative amplitude of ~ 2.7 percent and a pe- 2. Observations 
riod of ~ 6.4 days. 

We study the photospheric activity of CoRoT-6 by apply- CoRoT-6 has been observed during the second long run of 

ing the spot modelling approach already used for the high- CoRoT toward the 8 alactlc centre from 15 A P nl to 7 September 

precision light curves of CoRoT-2 (Lanza et al. 2009a), CoRoT- 2008 - Slnce the star ls relatively bright (V = 13.91 mag), the 

4 &za et al. 2009b), and CoRoT-7 (LanzaelalJ|2010). It al- tlme sampling is 32 s from the beginning of the observations. 

lows us to map the longitudinal distribution of photospheric ac- CoRoT performs aperture photometry with a fixed mask (see 

live regions and trace their evolution during the ~ 140 days of Fl 8- 3 ln [Fridlund etal. || 2010|) . The contaminating flux from 

CoRoT observations. The different rotation rates of active re- background stars falling inside the mask amounts to 2.8 ± 0.7 

percent and produces a dilution of the flux of CoRoT-6 that has 
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* Based on observations obtained with CoRoT, a space project op- been corrected by subtracting the 2.8 percent of its median flux 
erated by the French Space Agency, CNES, with partecipation of the value. The flux inside the star's mask is split along detector col- 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, umn boundaries into broad-band red, green, and blue channels. 
Germany, and Spain. To increase the signal-to-noise ratio and reduce systematic drifts 
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present in individual channels, we sum up the flux in the red, 
green, and blue channels to obtain a light curve in a spectral 
range extending from 300 to 1 100 nm. 

The obser v ations and data processing are described by 
iFridlund et al.l d2010l) . to whom we refer the reader for de- 
tails. The reduction pipeline applies corrections for the sky 
background and the pointing jitter of the satellite, which is 
particularly relevant during ingress and egress from the Earth 
shadow. Measurements during the crossing of the South Atlantic 
Anomaly of the Earth's magnetic field, which amounts to about 
7-9 percent of each satellite orbit, are discarded. High- 
frequency spikes due to cosmic ray hits and outliers are removed 
by applying a 5-point running mean. The final duty cycle of 
the 32-s observations is > 90 percent for the so-called N2 data 
time series that are accessible through the CoRoT Public Data 
Archive at IAS (http://idoc-corot.ias.u-psud.fr/). More informa- 
tion on th e instrument, i ts ope ration, and performance can be 
found in lAuvergne et al.l (120091) . 

The trans its are removed fro m the light curve using the 
ephemeris of IFridlund et al.l (1201 Ol) and the out-of-transit data 
are binned by computing average flux values along each orbital 
period of the satellite (6184 s). This has the advantage of remov- 
ing tiny systematic variations associated with the orbital motion 
of the satellite, still prese rving the informati o n on magnetic- 
induc ed variability (cf. lAlonso et al.l 120081: lAuvergne et alJ 
2009). In such a way, we obtain a light curve consisting of 1993 
mean points ranging from HJD 2454572.5075 to 2454717.4519, 
i.e., with a duration of 144.9443 days. The average standard er- 
ror of the mean points is 2.41 x 10~ 4 in relative flux units. We 
remove the long-term decrease of the flux, which may be due 
to the aging of the detector, by fitt ing a straight line to the data 
as in th e case of CoRoT-4 (cf. Aigrain et alj|2008c lLanza et all 
2009b). The maximum flux in the de-trended binned time se- 
ries at HJD 2454577.5245 is adopted as a reference unit level 
corresponding to the unspotted star, since the true value of the 
unspotted flux is unknown. 

3. Light curve modelling 

The reconstruction of the surface brightness distribution from 
the rotational modulation of the stellar flux is an ill-posed prob- 
lem, because the variation of the flux vs. rotational phase con- 
tains only information on the distribution of the brightness inho- 
mogeneities vs. longitude. The integration over the stellar disc 
effectively cancels any latitudinal information, particularly when 
the inclination of the rotation axis along the line of sight is close 
to 90°, as it is ass umed in the present case (see Sect. [4] and 
lLanza et al]|20( )9a). Therefore, we need to include a priori in- 
formation in the light curve inversion process to obtain a unique 
and stable map. This is done by computing a maximum entropy 
(hereafter ME) map, which has been proven to successfully re- 
produce active regio n distribution s and area variations in the case 
of the Sun (cf. lLanza et al1 l2007). Our approach has been com- 
pare d to other spot mod elling methods by iMosser et al.l (2009) 
and iHuber et al. I (120101) . The occultations of the stellar disc by 
the planet during transits can be used to resolve the fine structure 
of the spot pattern and study its evolution (e.g., Silva-Vali o et al.l 
2010). We plan to apply this approach to CoRoT-6 in a forth- 
coming paper (Silva-Valio et al., in prep.). 

In our model, the star is subdivided into 200 surface ele- 
ments, namely 200 squares of side 18°, with each element con- 
taining unperturbed photosphere, dark spots, and facular areas. 
The fraction of an element covered by dark spots is indicated by 
the filling factor /, the fractional area of the faculae is Qf, and 



the fractional area of the unperturbed photosphere is 1 — (Q+ 1)/. 
The contribution to the stellar flux coming from the k-th surface 
element at the time tj, where j = 1, ...,N is an index numbering 
the N points along the light curve, is given by: 



AF kj = / (wy){l-(e + l)/ + c 8 /+ 

2/[l + c f (l -fikj)])AkHkjW<jikj), 



(1) 



where /q is the specific intensity in the continuum of the unper- 
turbed photosphere at the isophotal wavelength of the observa- 
tions, c s and Cf a re the spot and facular contrasts, respectively (cf. 
Lanza et al. 2004), A k is the area of the fc-th surface element, 
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(2) 



is its visibility, and 



fiicj = cos xf/kj = sin i sin 9 k cos[4 + Q(fy - to)] + cos i cos 9k, (3) 

is the cosine of the angle iffy between the normal to the sur- 
face element and the direction of the observer, with i being the 
inclination of the stellar rotation axis along the line of sight, 
Ok the colatitude and ( k the longitude of the k-t\\ surface el- 
ement; Q denotes the angular velocity of rotation of the star 
(Q = 27r/P rot , with P tox the stellar rotation period), and to the ini- 
tial time. The specific intensity in the continuum v aries accord- 
ing to a quadratic limb-darkening law, as adopted bv lLanza et al.l 
(2003) for the case of the Sun, viz. /q oc « p + b p fi + c p fj. 2 . The 
stellar flux computed at the time tj is then: F (tj) = £a- &Fkj- To 
warrant a relative precision of the order of 10~ 5 in the compu- 
tation of the flux F, each surface element is further subdivided 
into 1° x l°-elements and their contributions, calculated accord- 
ing to Eq. ([TJ), are summed up at each given time to compute 
the contribution of the 18° x 18° surface element to which they 
belong. 

We fit the light curve by varying the value of / over the sur- 
face of the star, while Q is held constant. Even fixing the rotation 
period, the inclin ation, and the spot and facular contrasts (see 
lLanza et al]|2007l for details), the model has 200 free param- 
eters and suffers from non-uniqueness and instability. To find a 
unique and stable spot map, we apply ME regularization by min- 
imizing a functional Z, which is a linear combination of the x 1 
and the entropy functional S ; i.e., 



z= x 2 (f)-M(f), 



(4) 



where / is the vector of the filling factors of the surface ele- 
ments, A > a Lagrangian multiplier determining the trade- 
off between light curve fitting and regula rization; the expres- 
sion for S is given in lLanza et alJ (Il998l) . The entropy func- 
tional S attains its maximum value when the star is free of spots. 
Therefore, by increasing the Lagrangian multiplier A, we in- 
crease the weight of S in the model and the area of the spots is 
progressively reduced. This gives rise to systematically negative 
residuals between the observations and the best-fit model when 
A > 0. The optimal value of A depends on the information con- 
tent of the light curve, which in turn depends on the ratio of the 
amplitude of its rotational modulation to the average errors of its 
points. In the case of CoRoT-6, the amplitude of the rotational 
modulation is ~ 0.027, while the average errors of the points is 
~ 2.4 x 10" 4 in relative flux units, giving a signal-to-noise ratio 
of ~ 110. This is sufficient to adopt the sam e regularization cri - 
terion applied in the case of CoRoT-2 (see lLanza et al] |2009a). 
iterating the value of A until the condition |/i reg | =s e is met, 
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where /z reg is the mean of the residuals and eo - Co/ "V^V their 
standard error, i.e., the ratio of their standard deviation ctq in the 
case of the unregularized best fit (i.e., for X = 0) to the square 
root of the number N of points in each fitted subset of the light 
curve having a duration Aff (see below). 

In the case of the Sun, by assuming a fixed distribution of the 
filling factor, it is possible to obtain a good fit of the irradiance 
changes only for a limited time interval Aff, not exceeding 14 
days, which is the lifetime of the largest sunspot groups domi- 
nating the irradiance variation. In the case of other active stars, 
the value of Aff must be determined from the observations them- 
selves, looking for the longest time interval that allows a good 
fit with the applied model (see Sect.©. 

The optimal values of the spot and facular contrasts and of 
the facular-to-spotted area ratio Q in stellar active regions are 
unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product CfQ, so we can fix Cf and vary 
Q, estimating its best value by minimizing the x 2 of the model, 
as shown in Sect. [4] Since the number of free parameters of the 
ME model is large, for this specific application we make use of 
the model of Lanza et al. (2003), which fits the light curve by 
assuming only three active regions to model the rotational mod- 
ulation of the flux plus a uniformly distributed background to 
account for the variations of the mean light level. This procedure 
is the same used to fix the value of Q in the cases of CoRoT-2, 
CoRo T-4, and CoRoT-7 (cf. lLanzaeta l. 2009a b; L anzaetal.l 
l2010h . 

We assume an inclination of the rotation axis of CoRoT-6 
of i =s 89° in most of our models (see Sect. |4]l. Since the infor- 
mation on spot latitude that can be extracted from the rotational 
modulation of the flux for such a high inclination is negligible, 
the ME regularization virtually puts all the spots around the sub- 
observer latitude (i.e., 90° - i » 0°) to minimize their area and 
maximize the entropy. Therefore, we are limited to mapping only 
the distribution of the active regions vs. longitude, which can be 
done w ith a resolution better than ~ 50° (cf. [Lanza et al]|2007l 
2009a). Our ignorance of the true facular contribution may lead 
to systematic errors in the active region longitudes derived by 
our model because faculae produce an increase of the flux when 
they are close to the limb, leading to a systematic shift of the 
longitudes of the active regions used to reproduce t he observed 
flux modulation, as discussed by Lanza et alJ ( 20071) for the case 
of the Sun and illustrated bv lLanza et alJ (l2009al cf. Figs. 4 and 
5) for CoRoT-2. 



4. Model parameters 

The fundam e ntal stellar parameters are taken from 
iFridlund et ail (120101) and are listed in Table [Q The limb- 
darkening parameters q n , b p , and c p (cf. Sect. [3]) have been de- 
rived from Kurucz (2000) model atmospheres for r e ff = 6090 K, 
logg = 4.44 (cm s~ 2 ) and a metal abundance [M/H] = -0.2, 
by adopti ng the Co RoT white-band transmission profile given 
by lAuvergne et al.l d2009h . The limb-darkening parameters 
u + = 0.586 + 0.068 and u. = -0.12 + 0.13 obtained by fitting 
the transit cannot be used to constrain the above parameters 
because of their large errors (especially for uJ). However, they 
are compatible with those derived from Kurucz's models that 
yield u + = 0.646 and m_ = -0.110. 

The rotation period adopted for our spot modelling has been 
derived from a Lomb-Scargle periodogram analysis of the light 
curve giving P lot = 6.35 + 0.29 days (Scargle 1982). The uncer- 
tainty of the period comes from the total extension of the time 
series and represents an upper limit. As we shall see below, our 



Table 1. Parameters adopted for the light curve modelling of 
CoRoT-6. 



Parameter 



Ref." 



Star Mass (M ) 




1.05 


F10 


Star Radius (R Q ) 




1.025 


F10 


Teff (K) 






6090 


F10 


logg (cm s 


- 2 ) 




4.44 


F10 


Q p 






0.340 


LalO 


h 






1.024 


LalO 


Cp 






-0.378 


LalO 


P rot (days) 






6.35 


LalO 


e 




1.71 


xl0~ 4 


LalO 


Inclination 


(deg) 




89.07 


F10 


c s 






0.677 


LalO 


Cf 






0.115 


L04 


Q 






1.5 


LalO 


Af f (days) 




6.15 


, 7.247 


LalO 



" References: F10: IFridlund et all d2010h : L04: lLanza et alJ d2004h : 
LalO: present study. 



models show that the starspots have a remarkable differential ro- 
tation, which may contribute to an uncertainty of the stellar rota- 
tion period of w 10 percent, i.e., of ± 0.3 days, comparable with 
the above upper limit (cf. Sect. 15. 2\ . 

The polar flattening of the star due to the centrifugal poten- 
tial is computed in the Roche approximation with a rotation pe- 
riod of 6.35 days. The relative difference between the equatorial 
and the polar radii e is 1.71 x 10~ 4 which induces a completely 
negligible relative flux variation of » 10~ 6 for a spot coverage 
of ~ 2 percent as a consequence of the gravity darkening of the 
equatorial region of the star. 

The inclination of the stellar rotation axis has not been con- 
strained yet through the observation of the Rossiter-McLaughlin 
effect. Nevertheless, we assume that the stellar spin is normal to 
the orbital plane of the transiting planet, i.e., has an inclination of 
89°07 ± 0°30 from the line of sight, unless otherwise specified. 
This is compatible with the measured v sin i — 8.0 ± 1 .0 km s , 
the above rotation period, and a radius estimate of 1 .025 ± 0.026 
R H from stellar evol utionary models, that yield sin i = 0.98 (cf. 
IFridlund et al.l l2.010). Taking into account their errors, the mini- 
mum inclination still compatible with those stellar parameters is 
-50°. 

The maximum time interval that our model can accurately 
fit with a fixed distribution of active regions Aff has been deter- 
mined by dividing the total interval, T = 144.9443 days, into Nf 
equal segments, i.e., Aff = T/Nf, and looking for the minimum 
value of Nf that allows us a good fit of the light curve, as mea- 
sured by the x 1 statistics. We find that for Nf < 20 the quality 
of the best fit degrades significantly with respect to Nf > 20, 
owing to a substantial evolution of the pattern of surface bright- 
ness inhomogeneities. Therefore, we adopt Aff = 7.2472 days 
as the maximum time interval to be fitted with a fixed distri- 
bution of surface active regions in order to estimate the best 
value of the parameter Q (see below). For the spot modelling 
in Sects. IBTTI and 15.21 we shall consider the light curve interval 
from HJD 2454576.0 to 2454717.4519 because the light modu- 
lation is not well fitted during the first three days owing to the 
rapid evolution of some active regions (see Sect. 15.1b , and adopt 
N = 23 corresponding to Aff = 6.15 days, which provides a 
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better time resolution to study the evolution of the spot pattern 
during the intervals with faster modifications. 

To compute the spot contrast, we adopt the same mean tem- 
perature difference as derived for sunspot groups f rom their 
bolometric contrast, i.e., 560 K (Chap man et al.lfl994h . In other 
words, we assume a spot effective temperature of 5530 K, 
yielding a contrast c s — 0.677 in the CoRoT passband (cf. 
lLanza et al.l2007l) . A different spot contrast changes the absolute 
spot coverages, but does not significantly affect their longitudes 
and th eir time evolution, as discussed in detail by lLanza et al.l 
(2009j)- T he facular contrast is assumed to be solar-like with 
c f = 0.115 (lLanza etal.ll2004 . 

The best value of the area ratio Q of the faculae to the spots 
in each active region h as been estimated by means of the model 
of lLanzaetal] d2003l cf. Sect. [3]). In Fig. [Q we plot the ratio 
X 1 lx 2 ■ °f the total x 2 of the composite best fit of the entire time 
series to its minimum value v 1 . , versus Q, and indicate the 95 
percent confidence l evel as derived from the F-statistics (e.g., 
lLampton et~allll976l) . Choosing Af f = 7.2472 days, we fit the 
rotational modulation of the active regions for the longest time 
interval during which they remain stable, modelling both the flux 
increase due to the facular component when an active region is 
close to the limb as well as the flux decrease due to the dark 
spots when the same region transits across the central meridian 
of the disc. In such a way, a measure of the relative facular and 
spot contributions can be obtained leading to a fairly accurate 
estimate of Q. 

The best value of Q turns out to be Q a 1 .5, with an accept- 
able range extending from to ~ 4. The x 1 fluctuates owing to 
the numerical errors occurring during its minimization when a 
given best fit is computed. Therefore, it is its overall variation 
versus Q that is relevant here. In the case of the Sun, the best 
value of Q is 9 (lLanza et al.1 12007b . Thus our result indicates a 
lower relative contribution of the faculae to the light variation 
of CoRoT-6 than in the solar case. The relative amplitude of 
the rotational modulation of the star is ~ 0.027, showing that 
CoRoT-6 is more active than the Sun, for which a relative varia- 
tion of ~ 0.0035 is observed at the maximum of the 1 1-yr cycle. 
This may explain the reduced facular contribution t o the light 
variations of CoRoT-6, as suggested by Ra dick et al.l d!998l) and 
lLockwood et al.l (l2007h who find a smaller facular signature in 
stars more active than the Sun. 

The use of the chromatic information of the CoRoT light 
curves to estimate the spot and facular contrasts and filling fac- 
tors is made impossible by our ignorance of the unperturbed stel- 
lar flux levels in the different colour channels which are needed 
to disentangle the flux perturbations due to spots and faculae, 
respectively. The continuous variations of the observed fluxes 
do not allow us to fix such reference levels so that we cannot 
unambiguously attribute a given flux modulation to cool spots 
or bright faculae. Moreover, the oscillations of the fluxes in the 
individual colour channels on timescales comparable with the 
rotation period and the lack of precise passband profiles, make 
the extraction of information on spots and faculae unfeasible. 



5. Results 

5.1. Light curve model 

We initially apply the model of Sect.|3]to the whole out-of-transit 
CoRoT-6 light curve. Even without any regularization, the fit to 
the time interval between HJD 2454572.5075 and 2454576.0 
is not satisfactory owing to the remarkable flux variation on a 
timescale of ~ 2 days associated with the rapid evolution of 




6 
Q 



10 



12 



Fig. 1. The ratio of the x 2 of the composite best fit of the en- 
tire time series (from HJD 2454572.5075 to 2454717.4519) to 
its minimum value vs. the parameter Q, i.e., the ratio of the fac- 
ular area to the cool spot area in active regions. The horizontal 
dashed line indicates the 95 percent confidence level for^ 2 /^ 2 . , 
determining the interval of acceptable Q values. 
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Fig. 2. Upper panel: The out-of-transit light curve of CoRoT- 
6 (dots) in the time interval between HJD 2454572.5075 and 
2454582.0 and its unregularized best fit for a facular-to-spotted 
area ratio of Q - 1.5 (solid line). The flux is normalized to the 
maximum at HJD 2454577.5245. Lower panel: The correspond- 
ing residuals. Note the greater residuals before HJD 2454576.0 
(marked by the vertical dashed line) owing to the rapid evolution 
of some active regions. 



some active regions (cf. Fig. [2]). It cannot be adequately fitted by 
our model which assumes that the pattern of active regions stays 
constant during the time interval Aff > 2 days. Therefore, we de- 
cide to discard the first interval of the light curve, the endpoint 
of which is marked by the vertical dashed line in Fig. [2] and anal- 
yse the section from HJD 2454576.0 to 2454717.4519 covering 
141.4519days. We consider 23 individual intervals of A/f = 6.15 
days which warrant a good fitting and an adequate time resolu- 
tion to trace the evolution of the active regions. The best fit with- 



A. F. Lanza et al.: Activity, rotation, and star-planet interaction of CoRoT-6 



o 

CD 



1 -OOO <i 



0.995 



0.990 



0.985 



CO 
CD 



0.980 

0.004 

0.002 

O.OOO 

— 0.002 

-0.004 




4580 4600 4620 4640 4660 4680 4700 

Time (HJD — 2450000) 



Fig. 3. Upper panel: The out-of-transit light curve of C0R0T-6 (dots) and its ME-regularized best fit for a facular-to-spotted area 
ratio of Q — 1.5 (solid line) during the time interval between HJD 2454572.5075 and 2454717.4519. The flux is normalized to the 
maximum observed flux. Lower panel: The corresponding residuals. In both panels, the vertical dashed line marks HJD 2454576, 
i.e., the beginning of the interval considered for the modelling (see the text). 



out regularization {A = 0) has a mean p les - 5.739 x 10~ 6 and 
a standard deviation of the residuals <x = 4.155 x 10" 4 in rela- 
tive flux units which is 1 .7 times the standard error of the mean 
points. This is due to the flux fluctuations on a timescale of ~ 2 
days associated with a populations of short-lived active regions 
that can be only partially fitted by our approach (see below). The 
Lagrangian multiplier A is iteratively adjusted until the mean of 
the residuals p res - -4.962 x 1CT 5 - -cr / yffi, where N = 84 is 
the mean number of points in each fitted light curve interval Atf, 
The standard deviation of the residuals of the regularized best fit 
is a- = 4.728 x 1(T 4 . 

The composite best fit to the entire light curve is shown in the 
upper panel of Fig. [3] while the residuals are plotted in the lower 
panel. The fit is generally good, although the residuals some- 
times show significant excursions to be attributed to the rapidly 
evolving active regions with lifetimes of 1 - 3 days that cannot 
be fitted by our approach. In particular, during the first part of the 
observations up to HJD 2454576 marked by the vertical dashed 
line, the residuals are a factor of ~ 4 - 5 larger than during most 
of the subsequent interval. 

A ME model computed from HJD 2454576.0 to 
2454717.4519 with Q = 4 shows a marginally better fit 
with an average <x = 4.572 x 10~ 4 . Nevertheless, we decide to 
adopt Q — 1.5 in our analysis, corresponding to the minimum 
of the total x 1 of the fit obtained with the three-spot model, as 
described in Sect. [4] 

5.2. Longitude distribution of active regions and stellar 
differential rotation 

The distribution of the spot covering factor / versus longitude 
and time is plotted in Fig.|4]for our ME models with Q = 1.5. 
The origin of the longitude corresponds to the sub-observer point 
at the epoch HJD 2454576.0 and it increases in the same direc- 



tion of stellar rotation and the orbital motion of the planet. The 
lifetimes of the mapped active regions range from ~ 6 days, that 
is the time resolution of our spot modelling, to ~ 30 - 40 days, 
as can be inferred from the isocontours of their filling factors in 
Fig. [4] Most of the active regions grow and decay within five 
active longitudes that can be identified by their uniform migra- 
tion in the plot. Although there is always some degree of arbi- 
trareness in attributing a given active region to a specific active 
longitude, we trace the migration of these longitudes by white 
dashed lines. Their migration rates range from -3.88 + 0.75 
deg/day to 3 . 1 3 + 0.75 deg/day. If such different rates correspond 
to the angular velocity at different latitudes on a differentially 
rotating star, the relative amplitude of the differential rotation 
is estimated to be AQ/Q = 0.123 ± 0.020. This is actually a 
lower limit for the pole-equator angular velocity difference be- 
cause we do not have information on the latitude of the active 
regions which are probably confined to some range as we ob- 
serve in the Sun where spots appear only in the ±40° belt. A 
consequence of this interpretation is that the intersections be- 
tween the straight lines tracing the migration of different active 
longitudes in Fig.|4]do not necessarily imply a coincidence of the 
corresponding active regions because they are probably located 
at different latitudes. 



A ME spot model computed with Q = 4 leads to a somewhat 
smaller differential rotation, i.e., AO/Q =* 0.09 ± 0.02, because 
the reconstructed longitudes of the active regions are more af- 
fected by the flux modulation due to their facular component (cf. 
Sect. [3). Reducing the inclination i of the stellar rotation axis to 
70° does not significantly change the amplitude of the differen- 
tial rotation. On the other hand, for i = 50°, several small spots 
are ill-mapped making it difficult to trace the migration of the ac- 
tive longitudes rotating faster than the mean period of 6.35 days. 
This leads to an estimated AQ/O about half the above value. 
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Fig. 4. Isocontours of the ratio f/f mdx , where / is the spot cov- 
ering factor and / max = 0.0059 its maximum value, versus time 
and longitude for the ME models with Q - 1 .5. The two dashed 
black lines mark longitudes 0° and 360° beyond which the dis- 
tributions are repeated to easily follow spot migration. The con- 
tour levels are separated by A/ = 0.1/ max with yellow indicating 
the maximum covering factor and dark blue the minimum. The 
dashed white lines trace the migration of the active regions asso- 
ciated with each active longitude (see the text). 



The amplitude of the surface differential rotation of CoRoT- 
6 is in agreement with the mean value found for stars of com- 
parable effective temperature and rotation period by fitting the 
shape of their line profiles or by D oppler imaging techniques 
(iBarnes et al.l 120051: iReinersI 120061) . For CoRoT-4, that has a 
slightly hotter effective temperature of 6190 K and a mean ro- 
tation period of ~ 9 days, we estimate a lower limit for the sur- 
face differential rotation AQ /Q - 0.057 + 0.01 5 from the ME 
spot models with Q - 4.5 (lLanza et alJl2009a) . It may be ex- 
plained by assuming that the active regions mapped on that star 
are closer to the equator. By comparison, the relative angular ve- 
locity difference in the ±40° sunspot belt amounts to ~ 0.05. 

The total spotted area as derived from our modelling is plot- 
ted vs. time in Fig. [5] Since we removed the long-term linear 
trend in the flux, the associated variation of the spotted area is 
not included in our model. There is no clear periodicity in the 
variation, although timescales ranging from ~ 25 to ~ 60 days 
are apparent. 

5.3. Signatures of a possible star-planet interaction 

In the reference frame adopted for our spot modelling, which 
rotates with a period P rot = 6.35 days with respect to an iner- 
tial frame, the period between subsequent crossings of a given 
meridian by the subplanetary point is the synodic period P syn = 
22.25 days, defined as P~y n - P^ t —P~^, where the orbital period 
of the planet P 01 -b = 8.886 days. 

Considering the distribution of the spot filling factor / vs. 
time and longitude, we find a correlation of the longitude of 
some active regions with the meridian located at Al - -200° 
with respect to the subplanetary longitude (cf. Fig. [6). The white 
lines overplotted on the distribution of / show the motion of such 
a meridian vs. time with the period P syn . The active regions that 
reach a maximum filling factor close to their passage across this 
meridian are marked with crosses in Fig.|6]and listed in Tabled 
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Fig. 5. The total spotted area vs. time, as derived from our 
lightcurve model illustrated in Figs. [3] and [4] The uncertainty 
of the area (3cr errorbars) has been derived from the standard 
deviation of the residuals of the best fit to the light curve. 

where the columns from the left to the right report the time of 
their maximum filling factor, their longitude t u at that time, the 
filling factor per 18° longitude bin f/f max , normalized to the 
maximum filling factor over the whole map (/max = 0.0059), 
the distance 6€ of the active region from the active meridian, and 
the probability p of a chance association between them, respec- 
tively. The latter is computed according to: 



51 , 



(5) 



where n ar is the number of active regions appearing during a 
given synodic period and AL the longitude range where stel- 
lar activity is mainly localized. If activity were distributed uni- 
formly, AL = 360°, but in CoRoT-6 most of the active regions 
appear within a longitude range AL ~ 250°, owing to the pres- 
ence of active longitudes. Therefore, we adopt AL = 250° and 
War = 7 to compute the probabilities in the last column of Table|2] 
The adopted n m corresponds to the maximum number of active 
regions observed in one synodic period, thus yielding the maxi- 
mum probability of chance association. 

The total probability pj of a chance occurrence for the entire 
set of 13 active regions is: 



p T 



Up- 



min{p,-} 



(6) 



The division by minjp,} accounts for the additional degree of 
freedom related to the fact that the phase lag Al - -200° be- 
tween the active and the subplanetary meridians is not fixed a 
priori but is derived from the same set of data. In such a way, 
we obtain a probability of chance occurrence pj - 0.0075, indi- 
cating that the association of those active regions with the active 
meridian is significant. 

Spot models with a greater facular contribution, i.e., Q - 
4.0, confirm the presence of active regions associated with the 
orbital motion of the planet, although Al - -170°, owing to the 
systematic longitude shifts of the active regions associated with 
a greater value of Q. 

Since it is not possible to exclude a misalignment between 
the stellar spin and the orbital angular momentum, we have in- 
vestigated the impact of a different inclination of the stellar rota- 
tion axis on the above correlation finding that the active regions 
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Table 2. Stellar active regions associated with the active merid- 
ian. 
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listed in Table [2] are well reproduced also when i - 70°. On the 
other hand, only one third of them survive when i - 50°, which 
corresponds to the minimum inclination compatible with the es- 
timated stellar parameters (cf. Sect.©. In this case, no significant 
evidence of star-planet interaction is found. 

In Fig. UJ we plot the light curve of CoRoT-6 and mark the 
epochs corresponding to the passages of the active regions listed 
in Table|2]across the central meridian of the stellar disc. The two 
remarkable active regions reaching the maximum filling factor 
at HJD 2454583.68 and 2454710.53 correspond to flux minima 
at HJD 2454588.1 and 2454714.4, respectively, while the other 
smaller active regions are generally associated with changes of 
the slope of the light curve. Such slope changes can be attributed 
to the variation of their area (or filling factor) because the projec- 
tion factor is constant when an active region crosses the central 
meridian. 

The separations between the vertical marks in Fig.|7]depend 
on the epochs of maximum visibility of the corresponding ac- 
tive regions, which coincide with the times when they cross the 
central meridian of the disc and are in general not simultane- 
ous with the maxima of their filling factors as listed in Table [2] 
Specifically, given that the rotation period of our reference frame 
is P mi - 6.35 days, the mark corresponding to a given active 
region in Fig. [7] is shifted from the epoch of its maximum fill- 
ing factor by a time interval depending on its longitude that can 
range up to ±P m t. Their different longitudes and the fact that the 
active regions associated with the star-planet interaction do not 
appear periodically, lead to an uneven distribution of the corre- 
sponding marks vs. time in Fig. [7] 

6. Discussion 

CoRoT-6 is an interesting system to study possible interactions 
between the planet and its host star because its rotation is not 
synchronized with the planetary orbit allowing us to separate the 
rotational modulation from a signal at the orbital or the synodic 
period. For comparison, CoRoT-4 is a slightly hotter star, but it 
is synchronized, thus an activity signature of planetary origin is 
more difficult to find out. Nevertheless, a persistent active region 
close to the subplanetary longitude has been observed in that 



100 200 

Longitude (deg) 

Fig. 6. Isocontours of the spot filling factor as in Fig. [4] with 
white lines tracing the active meridian lagging the subplanetary 
meridian by M - -200°. White crosses indicate the relative 
maxima of the spot filling factor corresponding to the active re- 
gions listed in Table [2] 



star during the ~ 70 days of CoRoT observations (Lanza et al. 
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Fig. 7. The light curve of CoRoT-6 with the epochs of maximum 
visibility of the active regions listed in Table [2] marked by verti- 
cal dash-dotted lines. 
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l2009bl) . Another interesting case is that of r Bootis, also syn- 
chronized to the orbital motion of its hot Jupiter. In this star, 
Walk er et al. I ([2008) found evidence of a persistent photospheric 
active region leading the subplanetary longitude by « 70°. 

The other remarkably active star CoRoT-2 shows a short- 
term activity cycle with a period of ~ 28.9 days, as indicated by 
the variation of the total area of its spots, that is very close to 
ten synodic periods of the planet with respect to the mean stel- 
lar rotation period of 4.52 days. As a matter of fact, the synodic 
period of CoRoT-2 is only 2.84 days, which is shorter than the 
ti me resolution Aff = 3 .16 days of the spot mapping performed 
bv lLanzaet alJ (2009a) making it impossible to look for the ef- 
fects found in CoRoT-6 whose synodic period is 22.25 days, 
i.e., remarkably longer than th e time resolution of o ur mapping 
(Aff = 6.15 days). However. iPagano et al. I d2009l) find a sys- 
tematic increase of the variance of the optical flux immediately 
before the planetary transits in the light curve of CoRoT-2 which 
suggests short-lived active regions or flaring activity close to the 
subplanetary longitude. 

A major difficulty with the allegedly cases of magnetic star- 
planet interactions comes from the remarkably different be- 
haviours observed in different stars that makes it impossible to 
define a common and simple phenomenology. Moreover, in the 
case of CoRoT-6, there is no stable periodicity in the emergence 
of the active regions allegedly associated with the interaction 
and this makes the proposed correlation true only in a statisti- 
cal sense. The lifetimes of the active regions listed in Table|2]are 
comparable with the synodic period, thus, sooner or later, they 
should cross the straight lines marking the active meridian in 
Fig.|6]during their migration at the rates characteristic of the ac- 
tive longitudes to which they belong. Therefore, what suggests 
the presence of a star-planet interaction is not just the coinci- 
dence in longitude, but the fact that the maximum filling factors 
of those active regions are reached very close to such a crossing. 

This can be interpreted as the result of some triggering of 
activity associated with the passage of the active meridian over 
a previously existing active region, possibly due to an enhance- 
ment of the rate of magnetic flux emergence from the subpho- 
tospheric layers where magnetic flux is stored. In principle, two 
possible mechanisms can be invoked to explain this effect, re- 
spectively associated with the tidal force or with the magnetic 
field perturbation in the stellar corona. 

Tidal effects in CoRoT-6 are very weak because the semi- 
majo r axis of the planeta ry orbit is 17.9 times the radius of the 
star {Fridlund et al. 2010). Specifically, the tidal synchronization 
time for the star is ~ 900 Gyr for a tidal quality factor Q' - 10 6 , 
i.e., much long er than the estimate d stellar age ranging between 
1 and 3.3 Gyr (Frid lund et alJ 2010). Moreover, tidal effects are 
unlikely to be responsible for the observed star-planet interac- 
tion because they would produce two maxima of activity on op- 
posite hemispheres associated with the two opposite tidal bulges 
instead of the suggested intermittent activity associated with a 
meridian rotating with the synodic period. 

A perturbation of the hydromagnetic dynamo in the subsur- 
face lay ers of a star i nduced by a close-in planet has been consid- 
ered by Lanza (2008). The conjectured process assumes that the 
reconnection occurring in the stellar corona between the mag- 
netic fields of the star and the planet generates a perturbation of 
the magnetic helicity that is conducted down to the stellar sur- 
face and immediately below it along magnetic field lines. The 
longitude where the perturbation of the dynamo action is maxi- 
mum coincides with the footpoint of the magnetic field line con- 
necting the surface of the star to the planet. If the stellar dy- 
namo is not axisymmetric, when an active longitude crosses that 



footpoint the magnetic field in the subphotospheric layers can 
be intensified locally, giving rise to the observed maximum of 
activity in the active region. The phenomenon is not periodic 
because the emergence of the magnetic flux requires that the in- 
tensity of the unperturbed field be already close to the thresh- 
old for the buoyancy instability and this occurs only where there 
is a pre-existing active region. Since active regions appear in a 
non-periodic fashion and are preferentially located within active 
longitudes migrating at different rates, there is no well-defined 
periodicity in the star-planet interaction, but only some correla- 
tion with the activit y alrea dy found in the active longitudes. 

In the model of Lanza (2008), the coronal magnetic field has 
a non-potential configuration so its field lines are twisted and this 
can account for the obse rved lag Al be tween the active and the 
subplanetary meridians. lLanza I (12009) proposes that the mini- 
mum energy configuration of the stellar coronal field is a lin- 
ear force-free state, i.e., with the field satisfying the equation 
V x B — aB, where a is uniform. In principle, a unique linear 
force-free configuration can be found when the magnetic field 
components are given at the photosphere and the angle of twist 
of the field lines between the surface of the star and the orbit of 
the planet is specified. Unfortunately, in our case we have only 
information on the latter quantity, so there are potentially infi- 
nite configurations leading to the observed phase lag between 
the active longitude and the subplanetary point. For a phase lag 
Af = -200°, one of the possible configurations is obtained for 
a - 0.085, while the other field parameters specifying the varia- 
tion of the field strength with radius have been fixed at the values 
bo - —1.1 and cq = 1.0, as adopted by Lanza (2008). 

The topology of the coronal magnetic field plays an impor- 
tant role in the angular momentum loss of the star through its 
magnetized wind. The interaction between the planet and the 
coronal field could modify this topology leading to a reduction of 
the braking rate by * 30 percent (see Lanza 2010, for details). 
Indeed, the rotation period of CoRoT-6 seems to be quite short 
for a star of its mas s and a minimum estimated age of « 1 Gyr 
(Fridlund et al.ll2.Q10n. S pecifically, the gyrochronology relation- 
ship of Barnes (2009) gives a period of ~ 9 days, i.e., about 
30 percent longer than observed for a star of that mass and age. 
The difference is significant because the age spread allowed for 
^rot ~ 6 - 7 days is about ±100 Myr for a predicted age of ~ 500 
Myr. On the other hand, it is not possible to exclude that the 
star is remarkably younger than 1 Gyr because this lower limit 
was estimated from the Lithium abundance which has a weak 
dependence on the stellar age in that range of mass. Of course, 
a younger age (« 400 - 600 Myr) would reconcile gyrochronol- 
ogy with the measured rotation period. In the case of CoRoT-6, 
tidal effects on the evolution of the stellar spin are unlikely, given 
the very weak tidal force exerted by the planet that moves on a 
relatively wide orbit. 

Finally, we note that possible signatures of the star-planet in- 
teraction can be searched for also in the modulation o f the chro- 
mospheric line fluxes (e.g. Shkolnik et al. 2005, 2008). From the 
correlation amo ng chromospheric e mission, spectral type and ro- 
tation found by iNoves et aD (Il984l) . a residual flux in the core 
of the Ca II H&K lines R' HK » 2.8 that of the Sun can be es- 
timated (Foing, private comm.). It could be observable if high 
resolution (<: 60000) and sufficiently high signal-to-noise ratio 
(S/N Z 30 - 40) spectra of the Ca II H&K region are avail- 
able (cf., e.g., lReboloetal.lll989l) . Since the vsiru of CoRoT-6 
is 7.5 ±1.0 km s -1 , it could be possible to derive the longitudes 
of th e chromo s pheric plages by, e.g., the method proposed by 
IChar & Foing I (119931) . and study their possible relationship with 
the planetary orbital phase. 
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7. Conclusions 

We have applied a ME spot model to fit the optical light curve of 
the active planet-hosting star CoRoT-6 deriving maps of the dis- 
tributions of its active regions vs. longitude with a time resolu- 
tion of 6.15 days during 141 days of uninterrupted CoRoT obser- 
vations. Active regions are assumed to consist of cool spots and 
bright faculae with solar-like contrasts. The facular-to-spotted 
area ratio in each active region is assumed to be constant and 
the best value is found to be Q — 1 .5 with an acceptable range 
between and ~ 4. This is typical of a star with an activity level 
higher than the Sun. 

Our maps show five active longitudes where individual ac- 
tive regions preferentially grow and decay with typical lifetimes 
up to ~ 30 - 40 days. The migration of the active longitudes al- 
lows us to derive a lower limit for the surface differential rotation 
AQ./Q. = 0.12 ± 0.02. 

We find several active regions reaching a maximum of ac- 
tivity at a longitude lagging the subplanetary meridian by hi ~ 
-200°. A spurious coincidence, owing to a chance appearance 
of those active regions, can be excluded with a degree of con- 
fidence greater than 99 percent. This suggests the presence of 
some magnetic star-planet interaction, although the correlation 
we found is significant only in a statistical sense. 
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